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Abstract. We write down the anomalous dimensions of the fields of the 
Nambu-Jona-Lasinio model or chiral Gross Neveu model with a continuous 
global chiral symmetry for the two cases U{1) x U (1) and SU{M) x SU (M) 
at 0{1/N'^) in a expansion. 
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The Nambu-Jona-Lasinio, (NJL), model was introduced in Q as a the- 
ory with a continuous global chiral symmetry which is broken dynamically. 
It has also been studied in the context of hadronic physics following the early 
work of Q and lately in, for example, j^, and chiral perturbation theory Q 
since it corresponds to a low energy effective theory of the strong interac- 
tions. Further, there has also been significant interest in analysing models 
with four fermi interactions, like the NJL model, in other contexts such as 
the standard model, [||, |6|, ^. For instance, it was initially pointed out in |^] 
and subsequently studied in more detail in §, 0] that the Higgs boson of the 
standard model could be regarded as a composite field built out of fermions 
in much the same way that the a and vr fields of the two dimensional chiral 
Gross Neveu, (CGN), model are bound states of the fundamental fermions 
of the model, ||. 

One of the most widely used techniques to examine the CGN or NJL 
models is the large expansion where the number of fundamental fields 
N is allowed to become large, [^J, and can therefore be used as a di- 
mensionless coupling constant. Consequently one can show, for example, 
that the models are renormalizable in 2 < d < 4 dimensions in this ap- 
proach and simultaneously deduce that they possess a rich structure, 
such as dynamical symmetry breaking. Whilst the leading order large N 
analysis is well understood, it is important to go beyond the leading or- 
der to improve our knowledge of the quantum structure. Recently, new 
techniques to achieve this for theories with fermionic interactions were in- 
troduced in 1 10 1 for the 0{N) Gross Neveu model based on the critical point 
self consistency methods to calculate critical exponents in the bosonic 0{N) 
a model, [11, In this letter we present the results of the application of 



the same techniques to the 4-fermi models with continuous chiral symmetry 
by writing down the anomalous dimensions of each of the fields in arbitrary 
dimensions at 0{1/N'^). This is one order beyond any previous analysis and 
it is necessary to have such higher order expressions in order to improve our 
understanding of the areas mentioned above. A further motivation for such 
a computation lies in the accurate independent evaluation of these quan- 
tities in three dimensions in order to provide precise estimates to compare 
with recent numerical simulations of 4-fermi models, |jl3|. 

More concretely, the lagrangians of the models we have considered are, 

i, 

2r 
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for the NJL model with a global U{1) x U{1) chiral symmetry and 

for the same model but with an SU (M) x SU{M) chiral symmetry. The 
bosonic fields a and tt" are auxiliary and g is the perturbative coupling 
constant. In (2) the generalized Pauli matrices, A^j, 1 < a < (M^ — 1), 1 < 
I < M are normalized to Tr(A''A^) = 4T{R)6'''' and we used the conventions 



of |15| and [|^. Both (1) and (2) involve A^-tuplets of fermions 1 <i < N and 
this N will become our expansion parameter. We also use the conventions of 
[ll3| in defining the properties of 7^ as {7^^, 7^} = 0, tr(7^7'^i . . . 7^^") = and 
(7^)^ = 1 with trl = 2. We remark that our 7^ conventions in (i-dimensions 
retain the anticommutativity property. This differs from the definition given 



in 1 17] which is one formulation used to perform consistent renormalization 



calculations using dimensional regularization where the spacetime dimension 



is changed to provide a way of handling infinities, |17, 18|. There, |Tj 
one loses Lorentz invariance in the full d-dimensional space where [7^,7^] 
= when ^ is not an index in the physical spacetime. By contrast, the 
method of |0] uses propagators defined in arbitrary but fixed dimensions 
which are consistent with Lorentz and conformal symmetry. Therefore it 
seems more appropriate here to retain the anticommutativity of the 7^ for 



a Lorentz invariant formulation. (Indeed an alternative to [17| for treating 
7^ in dimensionally regularized calculations retains this condition, p^.) 
In the model with vr = was solved at 0(1/A^2) ^j^g gg^f 

consis- 
tency approach of obtaining critical exponents and we have used the same 
methods to treat (1) and (2). Briefly, this involves solving for critical expo- 
nents at the d-dimensional fixed point of the field theory where the model 
possesses an extra scaling or conformal symmetry. With this scaling prop- 
erty one writes down the most general form the propagators of the fields can 
take, consistent with Lorentz and conformal symmetry, where the powers of 
the scaling form are related to the dimension and therefore the anomalous 
dimension of the fields, |13]. To obtain analytic expressions for these critical 



exponents one represents the skeleton Dyson equations of various Green's 
functions by the scaling forms and solves the resulting representation of the 
Dyson equations for the unknown exponents order by order in The 
power of the method is illustrated by the fact that the leading order results 
are deduced algebraically and agree with previous work, whilst the new re- 
sults at the subsequent order are deduced with a minimal amount of effort, 
[10-12]. 
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For (1), the asymptotic scaling forms of the respective propagators are, 
in coordinate space, in the critical region, 

V'(^) - > ^(^) - J§y > - (3) 

where A, B and C are amplitudes which are independent of x and the 
dimensions of the fields are defined as 

a = ^ + 5^>/3=l-??-X<T,7 = l- ??-X7r (4) 

where d = 2^\s the spacetime dimension, r/ is the fermion anomalous dimen- 
sion and Xo- and Xtt are the anomalous dimensions of the respective vertices. 
The anomalous dimensions are expanded in powers of 1/iV via, for exam- 
ple, rj = X^i^i Vi/^^- By substituting the scaling forms (3) into the skeleton 
Dyson equations with dressed propagators, which are illustrated in figs. 1-3, 



we were able to determine r/2 using results obtained in |13|. We found 



m = 2^(2^ - 1) 

r(M-i)r(i-^)r(^ + i)r(^) ^> 



2 1 

^ A* + 7 + IT 

/i — 1 2/i 



(6) 



where ^'(^) = — 1) — "0(1) + "0(2 — A*) — i^ifJ-) and ^^{x) is the loga- 

rithmic derivative of the F-function, as well as Xo-i = Xni =0. Equation 
(6) represents the first 0(l/iV^) quantity to be determined for (1). 

The 0{1/N'^) corrections to the vertex anomalous dimensions were de- 
duced by following the analogous calculation for the 0{N) Gross Neveu 
model given in |jl6|. It involves studying the scaling properties of the 3- 
point functions at 0{1/N'^) using a method which developed the leading 
order work of ||T9| , ^ for the bosonic a model on CP{N). Rather than 
illustrate the large number of graphs which arise at 0{1/N'^), we have given 
in figure 4 the basic structure of the distinct graphs which arise, though the 
graphs with vertex counterterms have not been shown. The basic integrals 
corresponding to each of the graphs have been given in and it was there- 
fore a straightforward exercise to manipulate the graphs which occur in (1) 



to be proportional to integrals whose values are already known, |16]. For 
(1), at 0(l/iV^) we found that the degeneracy in the (cr, vr) sector was not 
lifted but unlike at leading order the fields now have a non-zero anomalous 
dimension, which is a new feature, ie 

X.2 - X.2 'li^^ 
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We have also carried out the same calculation for the non-abelian case 
(2) and we note the results obtained at leading order are 



^1 = TT 



m 

2 



1 , C2{R) 
M T{R) 



where fji = -2r(2^ - l)/[r(^ + l)r(^)r(l - ^)r(^ - 1)] and 



Xal 



XttI 



2(^ - 1) 



1 

M 



2(^ - 1) 

4C2{R)I, f'""^f^'^ : 



C2iR) 
T{R), 

C2{R) 1 
T{R) M 



C2(G) 



2T{R) 



(8) 

(9) 
(10) 



where A" A" 

C2{G) = M for SU{M), [|14|, |15|. Several leading order exponents were 



C2(G)5"' and C2{R) = (M^ - 1)/2M, 



calculated in [13| for SU{2) x SU{2) and our results are in agreement with 
them which provides us with a partial check on our calculation. At next to 
leading order the expressions are more involved compared to (6) and (7). 
We found 

2 



+ 




1 C2(R)\" ( , , 2 1 



1 C2(RW C2{G)C2{R) 

M T(R) ) 2T^{R) 



(11) 



2(/i - 1) 



Xa2 



and 



X7T2 



+ 



+ 



+ 



ml 



4(.-l)^ 

fiC2{R)C2{G) 
2T^{R) 

5fiC2iR) fl_ 
ip-l)TiR) [m 




T^{R) ; 



2(/i 



1 C2{R)y 

1) KM T{R) J 



^i{2i? -5^ + 4) 
(/i- 1)M 



1 3^2 (^) 
M T{R) 



(12) 



4(/. - 1) 



^(/i) + 



5 



+ 



+ 



+ 



(2^- 

I^C2{G) / C2{R) 
2T{R) \T{R) 

3^1^ (C2{R) 



1 , C2{R) 
M T{R) 



\T{R) 
2 C2{G) \ 
M 2T{R))_ 

1 C2{G) Y 
M 2T{R) ) 



1 C2iG) \ 
M 2T{R) ) 



4(^ - 1)3 V T{R) 



16(/x- 1)2M L V 
M2 



C2(G) 



+ T77 



1 

16 



1)3 



M J T{R) 
C2{G) Cl{G) 



+ 



3_ 
16 



(C2{R) 1 



T2(fl) 2T{R)M 8T'^{R) 
2 C2(G) 3C72(G)C2(i?) 
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2T^{R) 



+ 



2T^{R) 



M 



T{R) 
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AT{R) V T{R) 



. C2{G)\- 
2T{R)) 
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8 

\T{R) 
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+ 
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\T{R) 
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(ai-i)2; Vm2 

g2(g) \ ( C2{R) 
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+ 



G2{G) 

MT{R) ' MT{R) 
1 G2{G) Y 
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(13) 



where 0(/u) = '0'(/^) ~ V''(l)i ™d we have used the results of 10 in manip- 
ulating the f°'^'^ and d"'^'^ tensors which arise in the 3-loop graphs of fig. 4. 
We have expressed our results in as general a form as possible which allows 
one to check that each expression does agree with the analogous results of 
0, [l|] and the 0{N) model. 

We conclude with the observation that our results will prove to be ex- 
tremely useful in establishing which other models lie in the same universality 
class as (1) and (2) since, for example, we have provided independent ana- 
lytic expressions which can now be expanded in powers of e and compared 
with e-expansions of critical exponents of other models deduced from the 
corresponding perturbative renormalization group functions. 
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Figure Captions. 

Fig. 1. Dressed skeleton Dyson equation for ip. 
Fig. 2. Dressed skeleton Dyson equation for a. 
Fig. 3. Dressed skeleton Dyson equation for tt. 
Fig. 4. Vertex corrections with dressed propagators. 
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